Ground-state characterization of Nb charge-phase Josephson qubits 
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We present investigations of Josephson charge-phase qubits of SQUID-configuration inductively 
coupled to a radio-frequency-driven tank circuit enabling the readout of the states by measuring the 
Josephson inductance of the qubit. The circuits including junctions with linear dimensions of 60 nm 
x 60 nm and 80 nm x 80 nm are fabricated from Nb/A10 x /Nb trilayer and allow the determination 
of relevant sample parameters at liquid helium temperature. The observed partial suppression of the 
circulating supercurrent in the qubit loop at 4.2 K is explained within the framework of a quantum- 
statistical model. We have probed the ground-state properties of qubit structures with different 
ratios of the Josephson coupling to Coulomb charging energy at 20 mK, demonstrating both the 
magnetic control of phase and the electrostatic control of charge of the qubit. 
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I. INTRODUCTION 

Superconducting structures with mesoscopic Joseph- 
son tunnel junctions can provide a basis for electronic 
devices operating on single Cooper pairs. Prominent 
examples are the superconducting quantum bit (qubit) 
circuits which are regarded as promising elements for a 
scalable quantum computer T ne Josephson charge- 
phase qubit 0] is based on a Cooper-pair box Q of 
SQUID-configuration, i.e. a superconducting loop inter- 
rupted by two small-capacitance junctions with an island 
in between, which is capacitively coupled to a gate elec- 
trode (i. e. the Bloch transistor Q). The charging en- 
ergy Ec and the Josephson coupling energy Ej are typ- 
ically of the same order, so the dimensionless parameter 
A = Ej/Eq is of the order of one. Moreover, our cir- 
cuit comprises both the qubit and its readout [E]. The 
transistor can be operated as a box (qubit) whose dis- 
tinct quantum states with energies E n , n = 0, 1, 2, are 
associated with different Bloch-bands of the system 
The eigenf unctions \n, q) are the Bloch wave functions of 
a particle in the periodic (Josephson) potential. Here n 
is the band number and q the quasicharge governed by 
the gate voltage Vq, i. e. q = CqVg, where Cq is the 
gate capacitance, see, e.g., Ref. 0. The quantum states 
of the transistor also involve the phase coordinate <p set 
by the external magnetic flux <&d c applied to the SQUID 
loop. The variable <p behaves almost classically and is 
regarded as a parameter. Due to two control parame- 
ters (charge and phase), the eigenenergies E n (q,ip), the 
transition frequency uiq = (E\ — Eo)/h as well as - in 
the case of multiple qubits - the strength of mutual cou- 
pling, can be varied in the wide range. The read-out of 
our qubit can be performed similar as in the rf-SQUID- 
based impedance measurement technique pioneered by 
Rifkin and Deaver the qubit eigenstates can be dis- 
tinguished by the Josephson inductance of the effective 
weak link, i.e. the transistor, included in the loop whose 
impedance is probed by small rf-oscillations induced by 
an inductively coupled resonant tank circuit 0, , see 
the equivalent circuit in Fig. ^ 

Recently, Born et al. [T(| have confirmed the afore- 
mentioned read-out conception in the spectroscopic ex- 
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FIG. 1: Scheme of the experiment. The core element is 
a double Josephson junction (two crossed boxes) with a ca- 
pacitive gate coupled to its island, i.e. the Bloch transistor, 
embedded in a macroscopic superconducting loop. The loop 
is inductively coupled to an rf-driven tank circuit connected 
to a cold preamplifier. The capacitance between the loop and 
the ground is assumed to be much larger than the total ca- 
pacitance of the island. 



periment with Al qubit. The authors demonstrated the 
detection of microwave-power-induced interband transi- 
tions by tracking the resonant response of the tank. How- 
ever, due to the rather large ratio A = Ej/Eq ~ 30, 
manipulation of the qubit state via its gate was in their 
experiment limited to a narrow (about 5%) interval of 
values of the flux. In this paper, we present a compre- 
hensive mapping of the ground state of Nb qubits oper- 
ating in the charge-phase regime over a wide range of the 
parameters q and ip. This is done for the most illustra- 
tive case of the values of the parameter A ~ 1. In this 
regime, the behaviour of the system strongly depends on 
the value of A. Moreover, the effect of thermal fluctu- 
ations is studied comparatively at 20 mK and 4.2 K for 
the same samples. These investigations strongly benefit 
from the advanced technology for the fabrication of sub- 
orn Josephson junctions based on Nb/A10 x /Nb trilayers 
with a large superconducting energy gap A and a low 
sub-gap leakage current Applying this technology 

allows all inductive components of the circuitry shown in 
Fig. 2] to be integrated on one chip. 

The material presented in this paper is organized as fol- 
lows. In Sec. II, we briefly outline the considerations that 
motivated our choice of experimental parameters. In Sec. 
Ill, radio-frequency measurement results for niobium- 
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based Bloch transistors will be presented and quantita- 
tively analyzed. Especially the fluctuation-induced par- 
tial suppression of the critical current through Bloch 
transistors will be discussed within a quantum-statistical 
model. In Sec. IV, the concluding remarks will be made. 



II. DESIGN OF EXPERIMENT 
A. Parameters of the samples 

The operation of a charge-phase qubit bases on the 
single charge phenomena thus requiring the charging en- 
ergy Eq = e 2 /2C of the small island electrode (C = 
Cji + Cj2 + Cq denotes the corresponding total ca- 
pacitance including the capacitances of the junctions 
Cji,J2) to be large, compared to the thermal energy 
k-gT. In the case of our Nb technology the parameter 
Ec ~ 50 — 80 /j,eV roughly corresponds to temperature 
~ IK. The rather balanced situation with the target 
value A ~ 1, i.e. the Josephson coupling energy of in- 
dividual junctions Ejq — (^q/2tt)I c q ss 100 /ieV, where 
<I>o = h/2e w 2.07 ■ 10~ 15 Wb is the magnetic flux quan- 
tum, corresponds to nominal critical currents of the in- 
dividual junctions I c o of several tens of nA [l2| . 

The Bloch transistor is included in a double super- 
conducting loop of octagon shape with an outer dimen- 
sion of 1.0 mm (see Fig. 2). Such a gradiometer design 
improves the stability against homogeneous magnetic 
field noise. The total inductance value of such a loop 
L ps 0.7 nH was estimated using Mohan's formula |13j . 
This value is sufficiently small in the sense that, firstly, 
the corresponding magnetic energy is rather large, i.e. 
Em = $o/2L w 20meV and, secondly, the dimensionless 
screening parameter, i.e. (3l = 27rLi c o/$o <C 1, is small. 
The former relation ensures the suppression of flux fluc- 
tuations and fixes the phase tp across the transistor [l4|. 
The latter one ensures the required non-hysteretic opera- 
tion regime in the case of a single-junction SQUID circuit 
[1,^3 1 as we H as °f the realistic qubit circuit having finite 
asymmetry in the critical currents of individual junctions 

0. Moreover, the total inductance of the closed loop 
circuit is determined mostly by the transistor Josephson 
inductance. This inductance is directly related to the lo- 
cal curvature of the surface E n (q, if) taken for fixed q , 

1. e. 




(b) 



FIG. 2: (a) The image of the sample of the gradiometer con- 
figuration with the Bloch-transistor structure in the center, 
(b) The zoomed-in microphotograph of the transistor. 



The bare resonance frequency of such a tank circuit with 
quality factor Q w 250 is fo m 70 MHz. Due to the cou- 
pling to the qubit, the effective inductance of the circuit 
is changed |15| . 



(2) 



in a gate-charge and phase specific way. The resulting 
shift A/ of the resonant frequency, is [9j 
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Here, k = M/^/LLt < 1 denotes the coupling coefficient 
which is determined by the arrangement of the induc- 
tively coupled conductors. The integrated on-chip design 
of our qubit allows close mutual arrangement of the loop 
and the coil and, therefore, rather large values of k. The 
tank circuit is driven by a combined current signal, con- 
sisting of a dc part fixing the working point, and an rf 
component I r f at the frequency / close to fo- 

The cold preamplifier is based on an AGILENT ATF- 
10136 GaAs field-effect transistor. The output of the 
preamplifier is connected by the coaxial cable to the 
room-temperature amplifier. The overall gain of this two- 
stage scheme is about 26 dB in the range of typical reso- 
nant frequencies up to 90 MHz. The amplified rf voltage 
across the tank is fed to a lock-in amplifier measuring 
both the amplitude Vt and the phase shift a relative to 
the reference rf signal. 



In our measurements of the circuit in the ground state 
the band index n — 0, so we introduce the notation 
Lj(0,q, ip) = Lj(q, ip). 

The loop is coupled through the mutual inductance M 
to the coil of a resonant tank circuit formed by the dou- 
ble octagon-shaped spiral inductor with total inductance 
Lt = 0.15 — 0.2 //H (see Fig. 2a) and the capacitance of 
the coaxial cable connecting the tank to the preamplifier. 



B. Pre-characterization of Bloch transistors 

Embedding the Bloch transistors in a superconduct- 
ing loop, unfortunately, renders its pre-characterization 
in a simple dc configuration impossible. To get never- 
theless an estimate of the relevant parameters, i.e. the 
Josephson coupling energy Ejq and the charging energy 
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FIG. 3: (a) Typical dc I-V curve of the stand-alone test 
transistor from wafer B at 4.2 K. The superconducting gap- 
voltage in this sample is Vbap = 4.17 mV and the normal- 
state resistance ife = 67 kfi of the series double-junction 
(transistor), (b) Contour-plot of the differential conductance 
dl/dV in /iS of an isolated test transistor of wafer B as a 
function of gate voltage and bias voltage measured at 25 mK. 



Ec , we characterized the similar stand-alone Bloch tran- 
sistor structures fabricated from the same wafer, hence 
made from the same Nb/AlO^/Nb trilayer, and hav- 
ing the same dimensions. From I-V measurements at 
4.2 K (shown in Fig. we extract the sum asymp- 
totic resistance i?£ and the gap voltage Vcap- Resis- 
tance i?E is assumed to equal twice the normal state 
resistance i?N ~ 33.5 kfi of one from two nominally 
identical junctions while the superconductor energy gap 
A = eVcap/4 ~ 1.04 meV. Inserting these values into the 
Ambegaokar-Baratoff relation R^I c o ~ ( 7r /2)A/e (lt| . 
we estimate the critical current I c q and, finally, Ejq, the 
Josephson coupling energy of one junction as listed in 
Table I. Note that 7 c o values of 25 nA (wafer A) respec- 
tively 45 nA (wafer B) differ almost by a factor of 2. The 
radio-frequency measurements performed at 4.2 K with 
the single junctions inserted instead of the transistors in 
the similar loops (see Fig. [!}, i- e - the rf-SQUID configu- 
ration, gave almost similar values for I c q. 

The values of Ec were derived from the gate and bias 
voltage dependencies of the current of the stand-alone 
transistors at 20 mK and perpendicular magnetic field of 
up to 2T. This magnetic field partially suppresses the 
superconductivity of the Nb electrodes and enhances the 



single-electron tunneling at small voltage bias [17]]. For 
the bias below the gap voltage we found the characteris- 
tic diamond-like pattern resulting from Josephson quasi- 
particle cycles [18|; the bias voltage period provides an 
estimate of the charging energy (see Fig. and Table 
QJ , while the gate voltage period A Vq ~ 28 mV gives the 
value of the gate capacitance, Co = e/ AVc ~ 6 aF. 

TABLE I: Parameter overview for the experiments with the 
Bloch transistor samples Tl and T2. 

Tl (60 nm x 60 nm) T2 (80 nm x 80 nm) 



Ej (/LjeV) 


50 


95 


Ec (MeV) 


80 


45 


L (nH) 


0.7 


0.7 


Lt (mH) 


0.20 


0.15 


M (nH) 


4.6 


3.8 


k 


0.42 


0.39 


/c (nA) 


3.9 


17 


/fit 


0.008 


0.036 


fa (MHz) 


68.80 


76.18 


Q 


250 


235 



III. RADIO-FREQUENCY MEASUREMENTS 

A. Resonance curves 

First, the resonance measurements of the Bloch tran- 
sistors were performed in order to calibrate the rf oscil- 
lation amplitude of the phase, (p(t) = ifdc + a sin (27r/t), 
with a being proportional to the amplitude of the rf- 
flux oscillations in the loop, <& a = a$o/27r. The res- 
onance curves Vt — f for a moderate pumping ampli- 
tude are shown in Fig.^Ji. The curves exhibit character- 
istic "nodes" typical of rf-SQUIDs [l^, where, at certain 
pumping levels, the sensitivity of the tank oscillation am- 
plitude against the dc bias, i.e. the stationary Josephson 
phase difference tpd c , disappears. Assuming a sinusoidal 
current-phase relation (CPR), these "nodes" are related 
to the amplitude-dependent resonant frequency for fixed 
tpdc given by the formula |20| : 

/(a) = f ■ {1 + fc 2 /3 L [Ji(a)/a]co S¥ > dc } (4) 

with Ji(a) denoting the first order Bessel- function. If a 

equals one of the positive zeros cIq (i = 1,2,) of func- 
tion Ji, the resonant frequency is no longer dependent 
on ifdc or on the dc bias current. The amplitude for 
the first node can be assigned to the respective value a , 
thus providing calibration accor ding to a = 3.83V^/l4 1} , 
where V£ — maxVr- The good agreement between the 
calculated /(a)-dependence (dash-dotted line in Fig. 0Jl) 
and the position of the nodes illustrates, that the tran- 
sistor sample investigated behaves for rather large ampli- 
tudes approximately as an element with harmonic CPR. 
In other words, the contribution of higher harmonics in 
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between the driving signal J r f and the voltage oscillations 
is given by the formula 



;i=16nA,-\ <p c<c 71 
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FIG. 4: (a) Typical resonance curves for sample T2 at 20 
mK. Shown is the detected tank-circuit voltage as a function 
of the applied frequency of the rf-excitation for a driving- 
current of I T f « 320 nA. The dashed-dotted line shows the 
amplitude-dependency f(a) of the resonance-frequency for an 
assumed harmonic phase-dependence, as described by Eq. (|1J . 
(b) Typical resonance curves of sample T2 for a low excitation 
power, i.e. driving-currents I r { ~ 8 nA and I T { ~ 16 nA. 



the CPR to the resulting resonance frequency /(a) is 
small. 

For the small pumping levels, the resonance curve is 
shifted as a whole when changing the dc bias current, see 
Fig-Et- 111 contrast to the situation for the large ampli- 
tude of oscillations, the non-harmonic shape of the CPR 
plays a crucial role here as expressed by Eq. QJ relat- 
ing the local curvature of the band with the Josephson 
inductance. The frequency shift determined for the two 
opposite values of the external dc magnetic flux applied 
to the loop is no longer symmetric relative to the tank- 
circuit frequency since the reverse inductance LJ 1 (q,f) 
producing this shift (see Eq. (3)) is not exactly propor- 
tional to cos ip here. This asymmetric splitting appears 
so clear in our experiment because of the rather large 
value of the product k 2 Q(3h ~ 1.6 for sample T2 (cf. ex- 
periments on rf-SQUIDs with large values of this product 
in Ref. |2ll]). The minimum rf-excitation which allows us 
to work with a sufficiently high signal-to-noise ratio in 
our measurements corresponds to values a ss O.log • 

For a sufficiently small value of a the phase angle a 



tana = 2Q(A/// -0, 



(•5) 



with the resonance frequency detuning A/ resulting from 
the transistor's Josephson inductance Eq. and the rel- 
ative shift £ = (/ — /o)//o of the operation frequency /. 
The measurement of this angle for different values of the 
dc flux and gate charge at a sufficiently low temperature 
allows the curvature of the qubit ground state surface to 
be mapped. 



B. Radio-frequency measurements at 4.2 K 

As long as our samples are fabricated from Nb films 
having a critical temperature of about 9K, they preserve 
the Josephson properties and can be measured at a tem- 
perature of T = 4.2 K. Earlier, the rf measurements of 
transistors which comprise somewhat larger Nb junctions 
(with dimensions down to 300 run by 300 nm) and are in- 
cluded in low-inductance loo ps w ere successfully carried 
out at 4.2 K by IPichev et al. [22]. These samples showed 
a clear dependence of the phase angle on the applied dc 
flux and allowed the critical current of the junctions to 
be evaluated (about 55 nA). Due to the large capacitance 
of the island in these samples, the charging energy was 
small (about 4/xeV, i.e. 50 mK) and the values of pa- 
rameter A were rather large (> 30), so these transistors 
behaved classically at an operation temperature of 4.2 K, 
i.e. like two classical Josephson junctions connected in 
series. Surprisingly, our qubit samples with significant 
charging energy, i.e. A ~ 1, also exhibited a clear depen- 
dence of a on the dc flux at this temperature. 

Figure^ shows the periodic dependencies a-tpdc mea- 
sured in the qubit sample Tl at 4.2 K (solid line) and 
20 mK (dash-dotted line). For comparison, the corre- 
sponding curve (measured at 4.2 K) for the similar sin- 
gle Josephson junction included in the identical rf cir- 
cuit is shown by a dash-double-dotted line. This plot 
demonstrates first the reduction of the critical current 
of the transistor in the ground state at T = 20 mK 
(ksT <C Ej, Eq) in comparison to a single junction. This 
is due to the charging effect of the island (see, for exam- 
ple, the experiments |23l l24|). Secondly, one can see a 
further reduction of the critical current observed at the 
elevated temperature of 4.2 K. 

Such behavior can be explained by a simple model 
which takes the mixed Bloch states at values of thermal 
energy ksT into account which are comparable to the 
interband energies. In this case the observable supercur- 
rcnt (Is) is found as the quantum-statistical averaging 
over the canonical ensemble, 



h(<P,q) = (Is) = Tr(/ S e 



H/k-eT 



)/TrO 



,H/k B T 



), (6) 



where H denotes the Hamiltonian of the total system in- 
cluding the electromagnetic environment in the thermal 



5 



(a) 



E ~ 50 neV 



T1,T=20 mK- 
-T1, T=4.2 K 
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FIG. 5: (a) Phase-shift a as a function of the phase (ipdc oc 
&dc) for a single-junction (dash-double-dotted line) of wafer 
A at 4.2 K and for the double-junction sample Tl of wafer 
A at 20 mK (dash-dotted line) and 4.2 K (solid line). The 
nominal value of the Josephson energy of each junction of 
wafer A is about 50 /xeV. (b) Current-phase relations for the 
ground-band (dashed-dotted line), the first (dashed line) and 
second band (dotted line) of the Bloch-transistor calculated 
for E 30 = 100 fieV, E c = 45 ^eV and q = 0.5. The thick 
solid line indicates the observable value of the supercurrent 
7g (<,£>, g) derived from Eq. JHJ. 



equilibrium at a temperature T. In view of the small gate 
capacitance (Cq *C e 2 /2k^T) and the small loop induc- 
tance (L < ($o/27t) 2 /A;bT), the variables ip and q can be 
regarded as classical parameters. The full expression for 
the supercurrent operator Is was derived in Ref. Q . The 
diagonal matrix elements contributing to the expectation 
value entering in Eq. Jfjjl are, therefore, equal to 



(n\i s \n) 



2tt E n E J2 
$o Ej(<p) 



sirup {n\ cosx\ n ) 



(7) 



where the effective Josephson coupling energy of the tran- 
sistor is given by Ej((p) = (E^ + E 2 2 + 2Ej\Ej 2 cos ip) 1 ^ 2 
and x is the operator of the transistor island's phase (con- 
jugate to the operator of the island's charge). Note that 
the value | (n| cos xl 71 -) I < 1 ( see the plots of this matrix 
element for n — in Ref. 0) also depends on the phase 



ip, because ip determines the ratio A = Ej(ip)/ Ec- 

Note that due to the different signs of the term 
(n\ cosxjn) the supercurrent in the different Bloch bands 
(n\ Is\n), n — 0, 1, 2, ... has usually inverted the phase de- 
pendencies (see the calculated curves in Fig. |SJd as well 
as Fig. in Ref. [Bj]). These anti-phase terms contribute 
with the corresponding Boltzmann's factors to the ob- 
servable value Eq. ©. Finally, the supercurrent value to 
be detected in the 4.2 K measurement is equal to 

Is = (2Tr/<$>a)EnE J2 sm<p/Ej(<p) 

N 1 

(n\ cosx|n)e 



n— m— 
N 1 

/ EE< 

n— m— 



-E 11 (q-\-me,ip) 



(8) 



The summation over m takes into account both even and 
odd configurations of charge on the island. Both of these 
configurations are realized due to the unavoidable sin- 
gle electron tunneling at the elevated temperature. An 
analysis of Eq. JHJ shows that taking into account the 
three lowest energy bands, i.e. choosing N — 2, is suf- 
ficient for an adequate description of our experiments 
with the given sample parameters and temperatures up 
to T — 4.2 K. Then, higher bands (n > 2) are sparsely 
populated and do not essentially contribute to Is- As a 
tendency, the dependence of Is(<^> <z) on the gate charge q 
practically vanishes as soon as higher bands get involved. 
The phase dependence becomes almost harmonic. 

In Tabled we compare the critical-current values cal- 
culated according to our model, Ic th \ with those ex- 
tracted from the measurements, Ic , for both our tran- 
sistor samples at 20 mK and 4.2 K. From the experimen- 
tal data, we obtained the critical currents according to 
the relation 



(m) = $o tan(fo/2)q 
2tt fc 2 QLJi(a)' 



(9) 



following from Eqs. J2Jl and JSJ. Here we introduce the 
value 5a = max{ $Jc j a — min{$ dc j. a. For both temper- 
atures considered here we find that the theoretical and 
the experimental critical currents are in good agreement. 
Moreover, the fact that does not differ for zero tem- 
perature and 20 mK, indicates, that measurements at 
the base temperature of the dilution refrigerator explore 
indeed the ground state. 



C. Mapping of the ground state 

Figure |H1 presents the phase dependence as a function 
of both external dc-flux (corresponding to the Josephson 
phase ip) and the gate-voltage Vg proportional to the 
quasi-charge q for the samples Tl and T2. A nearly si- 
nusoidal dependence of the phase a on tp<ic is to be seen 
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TABLE II: Theoretical and experimental values of the critical 
currents of two Bloch transistor samples at a sufficiently low 
temperature of 20 mK and at 4.2 K. For the lji th ^ calculation 
based on the parameters Ejo and Eq listed in Table I, -Ejo = 
Eji = Ej2 is assumed as justified by the weak influence of 
the imbalance parameter b = (Ej 1 — Ej2)/(Eji + Ej2) found 
in the simulations. 
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FIG. 6: Phase-shift a as a function of the applied dc-fiux 
(oc Idc) and gate voltage Vg at T = 20 mK for samples (a) 
Tf and (b) T2. 



which is modulated periodically by the applied gate volt- 
age Vg- 

A closer look at Fig. reveals that the amplitude 
of the oscillations of a and, therefore, the critical cur- 
rent is smaller in Tl than in T2, as expected from the 
Ambcgaokar-Baratoff-values of Ej and also because of 
the smaller A-ratio, due to the charging effects. The nom- 
inal values Ejo differ by less than a factor of 2, whereas 
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FIG. 7: (a) Gate modulation-curves oi(Vg) of Tf for three 
different values of the applied magnetic flux. The curves are 
shifted for clarity, (b) Dependence of the peak-to-peak value 
Aa of the gate modulation as a function of the external mag- 
netic flux for Tl and T2. The negative values of Aa have the 
meaning of the reverse phase-dependencies in plot (a). 



the critical currents differ by more than a factor of 4. 
Related to that, the dc-bias dependence of a is closer to 
the sinusoidal one for sample Tl having the smaller value 
of A = 0.7 [H. 

Gate-modulation curves of sample Tl (low A-ratio) for 
different values ip^c are displayed in Fig.[7K- We find a 
periodic gate-modulation curve with a modulation depth 
|Aa| of 0.05 rad for (^dc = and of 0.09 rad for tpd c = w, 
whereas for <pdc = t/2, the gate charge sensitivity dis- 
appears almost completely. The gate-dependence period 
AVq is roughly 29 mV for both samples and complies 
with the value found for the stand-alone test transistor 
shown in Fig. |3Jd. The gate oscillations are le-periodic, 
since the periodicity does not change when applying a 
magnetic field of 2 T that is sufficient to cause intensive 
quasiparticle tunneling [171 ]. 

One should note that the gate modulation appears for 
sample Tl over the whole range of phase (fdc, hence al- 
lowing efficient qubit control over the whole flux bias 
range. On the other hand, for sample T2 (A = 1.9) the 
modulation depth is non-zero only in the vicinity of the 
point <^dc = 7T (cf. Ref. [l(i|[. where this range was no- 
tably smaller) and almost zero for the rest of the flux-bias 
range, see Fig. [Tjp. 
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FIG. 8: Experimental phase shift for samples (a) Tl and (b) 
T2 for two different values of the gate charge q — (open 
circles) and q = 0.5e (closed squares). The calculated depen- 
dencies are represented by solid lines (q = 0.5e) and dashed 
lines (q = 0), respectively. 

The surface plots of the phase shift presented in Fig.[f)| 
should reflect the dependence of the local curvature of 
the ground-states energy on q and <p (see Eq.^). How- 
ever, because of the finite amplitude of the oscillations, 

a ~ O.lag 1 ^, of phase tp these surfaces yield the values of 
curvature averaged over the finite interval. Such averag- 
ing of the reverse Josephson inductance is described by 
the integral 

1 f 1 dx 

LjHw)^- Lj l {q,<p + ax) (10) 

1" 7-i VI - x 1 

This expression makes it possible to compare the ob- 
tained experimental data with the corresponding depen- 
dencies following from the theory, taking into account 
the local curvature of the ground state energy and the fi- 
nite amplitude of the phase oscillations. By inserting the 
known parameters £jo, Ec, Q, A/, k and a into Eqs. Ipjjl. 
|[SJ and IjlUI) . we calculate the dc bias modulation curves 
for arbitrary charge on the transistor gate. The curves 
for q = and q = 0.5e are shown in Fig. |H1 These curves, 
which are based on input data partly deduced from the 
dc measurements described in Sec. II, agree well with 
the primary data from the rf experiments. 

IV. CONCLUSION 

The radio-frequency impedance measurements of the 
charge-phase qubit samples with balanced Josephson 



coupling to Coulomb charging energy ratio A clearly 
demonstrated the dependence of the curvature of the 
ground state energy on the control charge and phase in 
a wide range. The shape of the Josephson inductance 
surface of the transistor is well described by the Bloch 
band theory. 

An advantage of our Nb-technology is, that we were 
able to characterize our qubit system at a temperature 
of T = 4.2 K. In our investigation we found that the 
Josephson critical current of the tunnel junctions form- 
ing the qubit is scaled with their size, almost approaching 
the nominal Ambegaokar-Baratoff value. The availabil- 
ity of sub-0.1 /xm niobium-based Josephson junctions was 
extremely helpful as it offered the valuable new possibil- 
ity of studying the influence of thermal fluctuations in 
an extended temperature range without increasing the 
cryogenic efforts. The experimental results obtained at 
the large temperature can be interpreted within a simple 
quantum-statistical model of the Bloch transistor. 

Although the measured qubit samples Tl and T2 had 
a charging energy Ec of the island (equal to 80 /xeV and 
45 /ieV, respectively) much smaller than the value of the 
Nb energy gap, > 1 meV, this relation did not ensure the 
desired suppression of the quasiparticle tunneling [26| . 
As result, instead of the 2e-, the le-periodic dependence 
of the qubit Josephson inductance on the gate charge was 
observed. Such behaviour of the stand-alone Nb transis- 
tor samples was also observed in earlier measurements 
performed in dc configuration, see Ref. Illl and refer- 
ences therein. As was suggested in Ref. [llj, this be- 
havior is most probably due to possible intragap energy 
states formed in Nb/AlOx/Nb tunnel barriers or due to 
non-equilibrium quasi-particles in the outer electrodes of 
the transistors pa, |27|. Therefore, further improvement 
of the Nb technology of fabrication is required. On the 
other hand, the given Nb qubit samples can still operate 
in the "magic" points corresponding to the value of the 
control charge q = 0. The question as to the rate of the 
quasiparticle tunneling in the excited state, presenting 
the most critical mechanism of the qubit relaxation 
deserves a special study. 
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